A General Environment for the Treatment of Discrete Problems –

The slow convergence of the Incomplete Cholesky preconditioned Conjugate Gradient (CG) method, applied to solve the system representing a magnetostatic finite element model, is caused by the presence of a few little eigenvalues in the spectrum of the system matrix. The corresponding eigenvectors reflect large relative differences in permeability. A significant convergence improvement is achieved by supplying vectors that span approximately the partial eigenspace formed by the slowly converging eigenmodes, to a deflated version of the CG algorithm.